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Abstract. We constructed a local luminosity function (LF) of galaxies using a flux-limited sample (Si 70 > 0.195 Jy) 
of 55 galaxies at z < 0.3 taken from the ISO FIRBACK sur vey at 170 )im. The overall shape o f the 170-pm LF 
is found to be different from that of the total 60-fim LF iTakeuchi. Yoshikawa. fc Ishiil I2003T) : the bright end 
of the LF declines more steeply than that of the 60-pm LF. This behavior is quantitatively similar to the LF 
of the cool subsample of the IRAS PSCz galaxies. We also estimated the strength of the evolution of the LF 
by assuming the pure luminosity evolution (PLE): L(z) tx (1 + z)® . We obtained Q = 5.0^ 2 '5 which is similar 
to the value obtained by recent Spitzer observations, in spite of the limited sample size. Then, integrating over 
the 170-/xm LF, we obtained the local luminosity density at 170 /im, p_L(170/im). A direct integration of the 
LF gives pi(170/im) = 1.1 x 10 s h LoMpc -3 , whilst if we assume a strong PLE with Q = 5, the value is 
5.2 x 10 7 ft L Mpc" 3 . This is a considerable contribution to the local FIR luminosity density. By summing 
up with other available infrared data, we obtained the total dust luminosity density in the Local Universe, 
PL(dust) = 1.1 x 10 8 /i Z/0Mpc -3 . Using this value, we estimated the cosmic star formation rate (SFR) density 
hidden by dust in the Local Universe. We obtained psFR(dust) ~ 1.1— 1.2 h x 10 -2 Mq yr _1 Mpc -3 , which means 
that 58.5 % of the star formation is obscured by dust in the Local Universe. 

Key words, dust, extinction — galaxies: evolution — galaxy formation — galaxies: luminosity function, mass 
function — infrared: galaxies 



1. Introduction 

The luminosity function (LF) of galaxies is one of the 
fundamental statistics to describe the galaxy population 
in the universe. The far-infrared (FIR) LF is vitally im- 
portant to evaluate the amount of energy released via 
dust emission, and further, the frac tion of the star forma- 
tion activity hidden by dust (e.g.. IPerez-Gonzalez et alJ 
2005 [lLe Floc'h et alJl20o3 ITakeuchi. Buat. fc Burearellal 
2005). Not only the local LF but also its evolution plays 
a crucial role in understanding the cosmic star formation 
history. 

Most of the previous LF works at mid-infrared 
(MIR) and FIR wavelengths have been made based 
on IRAS da ta base 
Lawrence et alJ 
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The longest wavelength band of the IRAS is 100 /j,m. 
Subsequently, 15 a nd 90/xm LFs have been presente d 
based on ISO data llXull2000t rSerieant et allEoOll l2004h . 
Now, by the advent of Spitzer, 1 MIR(12 or 15 fim) LFs 
based on the 24-/zm band started to be available up t o 
z ~ 1 l|Perez-Gonzalez et alJEfiullLe Floc'h et alJl2005|) . 
Also recently, Spite er-based 60-p:m LF has been presented 
l|Fraver et al.ll2005|) . At 850 fxm, a LF of JRAS-selected 
sample of submill imeter galaxies has been published 
I Dunne et a"fll2000|h 

At wavelengths between 100 /jm and 850 /im, how- 
ever, only a very limited number of LFs have been stud- 
ied. Further, most of them are mad e from the sample 
select e d at shorter wavele ngths (e.g., iFranceschini et ail 
1998t iDunne et alJ l200fj|h or estimated/extrapolated 
from LFs of shorter w avelengths, e.g., 60 ^m (e.g., 
ISerieant fc Harrisonl I2005T) . A direct construction of the 
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LF is still rarely done up to now (see, lOvabu et al.ll2005h . 
Hence, it remains an important task to estimate the LF 
at wavelengths longer than 100 /im from a well-controlled 
deep survey sample. Wavelengths between 100 /im and 
850 /im are also very important in the context of the ex- 
tragalactic background radiation, particularly the cosmic 
infrared background (CIB). The CIB is now understood as 
an accumulation of radiation from dust in galaxies at var- 
ious redshifts (z). At the FIR, although the measured CIB 



is very strong (e.g., Gispert et alJl2000t lHauser fc Dwekl 
200ll iLaeache. Dole, fe PugetJl2005l among others), the 



properties of the sources contributing to the background 
is rather poorly known compared with other wavelengths. 
Thus, it is also of vital importance to have a LF at FIR 
for cosmological studies. 

In this work, we estimate the LF of the local galax- 
ies (z < 0.3) at 170 j tm based on the data obtained 
by FIRBACK survey l|Puget et alJ Il999j) . This paper 
is organized as follows: In Section [3 we describe the 
170 /im galaxy sample. We present the statistical esti- 
mation method of the LF in Section [21 In Section we 
show the LF and discuss its uncertainties. Section is 
devoted to our conclusions. We provide numerical tables 
of our LFs in Appendix A. Throughout this manuscript, 
we adopt a flat lambda-dominated cosmology with h = 
iVlOOfkms^Mpc- 1 ], and (Q ,A ) = (0.3,0.7), where 
J7o is the density parameter and Ao is the normalized 
cosmological constant. We denote the flux density at fre- 
quency v by S v , but for simplicity we use a symbol 5*170 to 
represent S v at a frequency (1.76 x 10 12 Hz) corresponding 
to A = 170 /Lim. 



2. Data 
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Fig. 1. The Kolmogorov-Smirnov test for the whole 
photometric sample and redshift subsample taken from 
FIRBACK 170 /im survey. The solid line shows the cu- 
mulative probability distribution of flux densities for the 
whole sample, while the dotted line depicts that of the 
redshift sample. The difference is found to be very small. 
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2.1. Parent sample 

The FIRBACK ( F ar-InfraRed BACKground) survey 
llPuget et all Il999t lLagache fc Doll l200lt iDole et all 

2001) is one of the deepest surveys performed at 170 /im 
by ISO using ISOPHOT llLemke et al.lll996h . It covers 
4 deg 2 on three fields. In this work, we use two of these 
fields: FIRBACK South Marano field (0.89 deg 2 ) and 
FIRBACK ELAIS Nl field (1.87 deg 2 ). 2 

The parent sample of FIRBACK is composed of the 
flux-limited sample of 141 sources with S170 > 135 mJy 
(3cr limit). Flux completeness of this parent sample is 
75 %, an d at flux density S170 ~ 200 mJy, it becomes 
- 90 % l|Dole et allboOlj) . We use the sources brighter 
than a flux density of 195 mJy in the following analysis. 



Fig. 2. The flux density-redshift distribution of our flux- 
limited sample. We used a subsample at z < 0.3 to con- 
struct the local luminosity function (LF). The vertical dot- 
ted line shows 6*170 = 0.195 Jy, and the horizontal dashed 
line depicts z — 0.3. 

2.2. Redshifts and completeness 

Redshifts are measured for 58 galaxies out of 69 galaxies 
above the flux density of 195 mJy, i.e., the redshift com- 

2 The choice of the two fields is due to the fol low-up allotment 
was different for ELAIS Nl /South Marano iDennefeld et alJ 
l2005ft and for ELAIS N2 jTavlor et al.ll2005l) in the FIRBACK 
project. Consequently the conditions of the data aquisition are 
different between them. A coherent treatment remains as a 
future work. 
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pleteness of the sample used is 84 %. The redshift measure - 
ments have been perfo rmed bv | (^ hanman et alJ l|2002h , 
IPatris et al] l|2003|) . and lDennefeld et all l)2005|) . Since the 
redshift measurement becomes more difficult toward the 
fainter sources, the completeness depends systematically 
on the flux levels. Thus, we should examine whether the 
redshift selection distort the flux distribution of the sam- 
ple. 

We performed the Kolmogorov -Smirnov test (e.g., 
lHoelll97lllHaiek. Sidak. fc Senll999(l to compare the flux- 
limited sample (Siro > 195 mJy) with the redshift sample 
(see Figure The maximum difference between the cu- 
mulative distribution functions of flux-limited and redshift 
samples are 0.0377. This shows that we cannot reject the 
null hypothesis that the two samples are taken from the 
same parent distribution. Hence, we use the redshift sam- 
ple as an unbiased subsample of the whole flux-limited 
sample and simply multiply the inverse of the complete- 
ness to obtain the final galaxy density. The distribution of 
the flux densities and redshifts of the sample is shown in 
Figure 

The redshift comple teness of the sample is also tested 
by V /Vma* statistics l|Schmidtl Il968t iRowan-Robinsonl 
Il968|) . Here V is the volume enclosed in a sphere whose ra- 
dius is the distance of a considered source, and Vmax is the 
volume enclosed in a sphere whose radius is the largest dis- 
tance at which the source can be detected. If the sample is 
complete, V/V max values of the sample galaxies is expected 
to distribute uniformely between and 1, with an average 



-1/2 



n: 



(V/Vmax) — 0.5 and a standard deviation (12n) 
sample size). For our redshift sample, the mean and stan- 
dard deviation of the V/V maj ^ is 0.66±0.23, i.e., the sample 
can be regarded as complete. Moreover, it is larger than 
0.5 (but within the un certainty) , sug gesting; the existence 
of evolution (see, e.g.. |PeacockHl999l p. 444). 

For the estimation of the local LF, we use a subsample 
of galaxies with z < 0.3. The size of this 'low-z' subsample 
is 55. The mean and median redshift of this low-z sample 
is 0.12 and 0.09, respectively. 

3. Analysis 

3.1. K -correction 

The monochromatic luminosity at observed frequency f bs 
is obtained by 

47rd L (z) 2 S„ ob3 



-(l+z)i/ ob3 



1 + 2 



(1) 



where C v is an energy emitted per unit time at frequency 
v, di,(z) is the luminosity distance corresponding to a red- 
shift z, and ^ bs and v cm are observed and emitted fre- 
quencies, respectively. In order to estimate the luminosity 
function at 170 fim, the if-correction is required. 

However, the amount of the i^-correction may be un- 
certain, because the present sample is observed at one 
waveband. If we assume a 'cool' dust galaxy, the spec- 
tral energy distribution (SED) rises toward longer wave- 
lengths, whilst it decreases if we adopt a starburst SED 



faggfi.g..lTak'eiic}ii et al.ll200llhllLa g ache. Dole, fc PngeJ 
2003). To explore the effect of the fsT-correction, we use a 
power-law approximation with the form of 



(2) 



As for 3, we consider 3 — 1.0 (starburst galaxies), 0.0 (in- 
termediate galaxies), and —0.5 (cool galaxies). We adopt 
these values according t o the phenomenological l y con- 
structed model SEDs of iLaeache. Dole, fc Pugetl l|2003h 
(see their Figure 4) . Then, the luminosity at the observed 
frequency ^ b s becomes 



^fobs — ^Obs-^obs 

L^Jl + z)-^ 
4TTd L (z) 2 iy ohs S^ bs (l + z)-^ 



(3) 



By the same manner, the limiting luminosity of a survey 
with flux density detection limit S^ m depends on the SED 
via 3, 



Lj™ = 47rd L (z)^ obs 5^ a (l + z)-W+V . 



(4) 



This is shown with the present sample in Figure The 
luminosity L\jq is that at emitted wavelength of 170 /an, 
i.e., L cm measured at 170 ^m. As we see in the followings, 
this dependence of the limiting luminosity on 3 potentially 
affects the estimation of the LF. 

3.2. Estimation of the luminosity function 

We define the luminosity function as a number density 
of galaxies whose luminosity lies between a logarithmic 
interval [log L, log L + d log L] : 3 



4>{L) 



dn 



d log L 



(5) 



In this work, we denote the luminosity at a certain fre- 
quency or wavelength as L u = vC v . 

3.2.1. Parametric estimation 

First we performed a parametric maximum likelihood es - 
timation of the LF l|Sandage. Tammann. fc 
Note that we can do this analysis directly on the data, 
being independent of the nonparametric result, i.e., this 
is not a fitting to the nonparametric LF. It is known that 
the 60-Atm LF i s wel l expressed by a function given by 
Sau nders et al which is defined as 



exp 



1 

2^ 



L 



(G) 



where log 2 x = (log a:) 2 . Since various LFs can be approx- 
imated by this functional form, we adopt Equation (jSJ in 
this work. 



We denote log x = log 10 x and In x 



log e x, respectively. 
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Fig. 3. The luminosity-redshift distribution of our redshift sample taken from FIRBACK 170 /im survey. The lumi- 
nosity is that at emitted wavelength of 170/mi. The dotted curves represent the limiting luminosities corresponding to 
the flux density detection limit of 195 mJy at each redshift, including the effect of the if-correction. The power-law 
index of the spectral energy distribution (SED), (3 is 1.0, 0.0, and —0.5 from left to right (see the main text). Vertical 
dashed lines show z = 0.3, which is used to define our low-z sample. 



We use Equation © for the parametric maximum like- 
lihood estimation. Then the likelihood M. is expressed as 



JV 

n 



M{L <fl a 1 (j\{L l ,z i } l 



)=n 



log L„ 



<j}{L)d\ogL 



— exp 



log L„ 



exp 



1 

2^2 



log 2 1 



L 



Note that, in principle, the parametric estimation proce- 
dure is dependent on (3. We can obtain the parameters of 
the LF by maximizing Equation l|7|l with respect to a, 
and a llSandaee. Tammann. k, Yahijl979llSaunders et alJ 
ITflOollTakeuchi. Yoshikawa. fc Tshiill200.^. 



where 



t . . — r lim (vA 



dlogL 

However, because of the small size of the present sam- 
pie 1 and relatively narrow range of their luminosity, it is 
difficult to put a reasonable constraint to the faint-end 
slope of the LF. Hence instead, as we explain later, we 
(8) assume a certain value for the faint-end slope a. 
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Fig. 4. The logarithmic likelihood hi Ad for the parameter estimation of L* [Lq] and a when we fix a — 1.25. The 
outermost contours indicate the 68 % confidence level. 



3.2.2. Nonparametric estimation 

We also estimate the LF nonparametric ally via an 
impro ved version of the C~ method of iLvnden-Belll 
I l97ll). implemented to have the density normalization 
( Choloniewskil Il987|) . This method is a kind of maxi- 
mum likelihood methods insensitive to the density fluctu- 
ation. This method and its extension are fully described 
and c arefully examined by iTakeuchi. Yoshikawa. fc Ishiil 
ll2000lh 4 

We note that the SED slope (3 also affects the 
nonparametric estimation of the LF. In the case of 
C~ method, the definition of C~ includes L 1 ) ™ (see 
Figure 2 of ITakeuchi. Yoshikawa. fc Ishiil l200(TK Thus, 
it will be important to explore the systematic effect in- 
troduced by if-correction. The uncertainty (68 % con- 
fidence limit) is estimated by the bootstrap resampling 
ijTakeuchi. Yoshikawa. fc Ishiil2 000). Additionally, we also 
estimated the LF and uncertainty including the observa- 
tional measurement errors. For the density normalization, 
we took i nto account the source extraction completeness 
(~ 90 %: iDole et al~ll200l[) and the redshift measurement 
(84 %). 

4. Results 

4.1. Parametric result 

We fixed the faint-end slope of the FIRBACK 170-/im LF 
to be 1.25. This is very close to that of the total 60-/im 
LF of the IRAS PSCz sample, and the same as that o f its 
'cool' subsample l|Takeuchi. Yoshikawa. fc Ishiill2003lh 

4 We found that the other density-insensitive nonpara- 
metric estimators discussed in ITakeuchi. Yoshikawa. fc Ishii 
( 2000) were not very suitable for the present small sa mple 
analysis: Both of the methods of ICholoniewskil il986t) and 
lEfstathiou et, all (|l98ST) need to divide the sample into small 
bins. For the present sample (55 galaxies), we could not find 
stable solutions for these estimators. 



We obtained L« = 1.1 x 10 9 /i~ 2 [£©] (with a 68-% 
confidence range of 3.0 x 10 8 /i~ 2 -2.3 x 10 9 h~ 2 [L©]) and 
g = 0.41 (with a 68-% confidence range of 0.35-0.50). 
The likelihood contours are shown in Figure The outer- 
most contours indicate A In Ad = hiAd— lnAd max = —0.5, 
corresponding to the 68-% confidence limit. In Figure 0] 
we present In Ad defined by Equation 0, as a function 
of (L*,er). These two parameters are rather strongly de- 
pendent with one another, and as a result, the contour 
is elongated along with the diagonal direction in each 
panel. The density normalization was 4>* — (1-0 ± 0.4) x 
lO^/i 3 [Mpc" 3 ]. As seen in Figure 0] the result is almost 
independent of the assumed (5. The result is also found to 
be quite robust against the value of a in a plausible range 
of a= 1.1-1.3. 

ITakeuchi. Yoshikawa. k. Ishii ( 2003) presented the pa- 
rameters for the LF of the IRAS PSCz sample. The pa- 
rameters for the LF of the whole sample are (a, L*,cr) = 
(1.23,4.34 x 10 s h- 2 [L Q ], 0.724). Clearly, the parame- 
ters for the 170-/zm sample are different from these val- 
ues. Particularly, a which determines the steepness of the 
bright end is significantly smaller than the total IRAS LF. 
From the above likelihood analysis, a — 0.724 was rejected 
with a confidence level of more than 99.9 %, even for the 
small number of galaxies. That is, the bright end of the 
present LF declines more steeply than that of the total 
IRAS LF. 



4.2. Nonparametric result 

Nonparametric LF estimates are presented in Figure |5| 
Symbols are the LFs obtained by the C~-method. Vertical 
error bars show the 68 % uncertainty, obtained by boot- 
strap resampling of 10 times. Open squares, open trian- 
gles, and open diamonds represent (5 = 1.0, 0.0, and —0.5, 
respectively. Solid lines are the analytic expression of the 
LF with the parameters we obtained in Section 14.11 We 
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Fig. 5. The local luminosity function of ISO FIRBACK 170 /xm galaxy sample. Solid curves show the parametric 
form estimated from Equation J7J. Symbols represent the C~ nonparametric LFs, respectively. Open squares, open 
triangles, and open diamonds are the LFs adopting (3 = 1.0, 0.0, and —0.5, respectively. Vertical error bars are 68-% 
confidence ranges. Left panel shows the LF from the original data, while the right panel shows the LF from the data 
convolved with observational measurement errors. In the left panel, horizontal bars simply represent the bin width 
(0.3 dex). In the right panel, in contrast, horizontal bars are the convolution of the bin width and the luminosity 
uncertainty introduced by the photometric error. For visual simplicity, we put horizontal bars only on the case of 
j3 = 0.0 in each panel. 



present the numerical tables of our C LF estimates also 
in Appendix A. 

The left panel is the LF estimates from the original 
data themselves, while the right panel is the ones from the 
data convolved with measurement errors. We performed 
the error convolution procedure assuming a Gaussian dis- 
tribution for the measurement errors in flux density, with 
quoted flux density errors as standard deviations for each 
data. This procedure slightly blurs the LFs horizontally. 
As a result, we have an additional bin at the highest lu- 
minosity, and the 68 % uncertainty levels are broadened. 
However, the effect is rather small, and the LF estimates 
are quite robust except for the highest-luminosity unstable 
bins. 

We see that different values for (3 do not affect the re- 
sult very strongly, but at the highest-luminosity bins, sys- 
tematic effects are relatively large, a factor of 0.5-0.7 dex. 
In contrast to the large effect of (3 we found in Figure [3] 
the LF estimates are quite robust against (3. Hence, here- 
after, we do not show all the results with respect to (3, but 
only restrict ourselves to the results with (3 = 0.0 without 
loss of generality. 

We should note the upward deviation of the C~ LF 
from the analytic result. Although the error bar is large, 
the C~ estimates (open symbols in Figure EJ) are about 
an order of magnitude larger than that of the analytic 
value. It is worth examining if this deviation is real or 
merely an artifact of the poor statistics. To see this, we 



show the luminosity distribution of our present sample in 
Figure Recall that this sample consists of galaxies only 
at z < 0.3. We see an excess at the highest luminosity bin. 
Examining the subsamples, we also find a similar excess 
in the ELAIS Nl field. We, however, also see that the 
sample in South Marano field is interesting: the galaxies 
in this subsample are biased toward higher luminosities. 
The superposition of these effects makes the brightest end 
of the nonparametric LF deviated from the analytic one. 

Why does the analytic maximum likelihood not reflect 
this sample property? We see that, though there is an 
excess at L170 ~ 10 n h Lq, the majority of the sample is 
located around 10 10 h L© (the peak in Figure EJ. Hence, 
the parameter estimation is practically controlled by the 
sample around L*, and the bright galaxies could hardly 
give a strong effect to the final estimation. In addition, 
the peak is dominated by the sample from ELAIS Nl field, 
and the peculiar luminosity distribution of South Marano 
field affected the result little. We, hence, should not overly 
rely on the analytic result and parameters, but rather we 
should use the nonparametric C~ LF directly. These are 
small sample effects, and we should wait for the larger 
sample to clarify the detailed shape of the LF. 

We show the redshift distribution of the present sample 
in Figure [7| The histograms show the distribution of the 
present data. The thick solid curve is the expected redshift 
distribution of galaxies calculated from the nonparametric 
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Fig. 7. The redshift distribution of our flux-limited sample. The dotted histogram is the distribution of the whole 
sample, the dashed histogram shows the the subsample in South Marano field, and dot-dashed one presents the 
subsample in ELAIS Nl field. These histograms present the number of galaxies in each bin (Az = 0.04). The thick 
solid curve is the expected distribution of galaxies calculated from the nonparametric local LF (Figure EJ. The thin 
solid curve is the expected distribution for the pure luminosity evolution with L(z) oc (1 + z) 5 . 



local LF we have obtained. This is calculated as follows 



^ — ^ L survey 



logL min (z) 



d 2 V 

( t > { L ™)-T7Ti-d\ogL in , (9) 
di Idz 



where L m - m (z) is the minimum detectable luminosity 
Ll m (z) [see, Eq. JHJ], and 



d 2 V 



d-L(z) 



dfldz H (l + z)V^o(l + 2) 3 + A (10) 

for th e flat lambda-dominated universe (see, e.g.. |Peeblesl 
Il993lh i.e., {d 2 V/dridz)dz is the comoving volume between 
[z, z + dz] per unit solid angle. This curve is calculated to 
apply to the solid angle of the survey area by multiplying 
^survey = 8.4 x 10~ 4 sr. We used the nonparametric LF 
for Figure [7| because the exact shape of the bright end of 
the LF is crucial to examine the tail of the redshift distri- 
bution of the source toward higher z, and as already dis- 
cussed, the analytic function underproduces galaxies with 
Lno ^ 3 x 10 10 Lq. We find a density excess at z ~ 0.2 in 
Figure[7| 5 However, apart from the bump, both fields have 
long tails toward higher redshifts than expected from the 



5 If we use the classical 1/ Vmax-method l)Schmidtll 968: Ealcs 
Il993f) . the estimator is affected by this bump and results in 
a (fake) flatter bright end of t he LF. This is a well-known 
draw back of the 1/Vmax-method ijTakeuchi. Yoshikawa. fc Ishii 
2000;, and references therein), and we must be very careful 
about its usage. We address this problem in Appendix B. 



local LF. It may suggest the existence of galaxy evolution. 
We explore the effect of the evolution in Section l5~2*l 



5. Discussion 

5.1. Shape of the 170 \xm LF 

The overall shape of the 170-^m LF of the FIRBACK 
galaxies is different from the IRAS 60-^m LF. As we 
discussed in Section 14.11 the parameters of the analytic 
LF suggest a steeper slope for the bright end of the LF. 
Although the nonparametric LF revealed that the bright- 
est part of the LF is significantly higher than that of the 
analytic LF, the 170-/xm LF decreases more rapidly from 
the knee of the LF to the bright end than the IRAS 60-/im 
LF. Here we compare the 170-/mi LF with other FIR LFs 

ob tained to date. 

iTakeuchi. Yoshikawa. fc Ishiil <|2003() divided the IRAS 
sample into two categories, warm and cool subsamples, 
using the flux density ratio criterion of <Sioo/<S'60 = 2.1. 
The parameters of the cool galaxies are (a,L*,cr) = 
(1.25,9.55 x 10 9 h- 2 [i ],O.5O). The a is almost the same 
as that of the total IRAS LF within the quoted error, but 
L* and a are much closer to those of our present LF. The 
resemblance of our LF to the LF of cool IRAS galaxies 
may be reasonable, since the present sample is selected at 
170 /im band, where we can detect cooler galaxies more ef- 
fectively than 60 /im. ISoifer k. Neueebauerl l)l99l|) derived 



8 



T. T. Takeuchi et al.: 170 /im galaxy luminosity function 



25 [ 



20 - 



ISO FIRBACK 170/xm 
(|S=0.0) 



S 15 



tt 10- 

XI 

a 

3 

5 - 



Whole sample 
South Marano 
ELAIS Nl 



OIL 
10 9 



10 10 10" 
\70fj,m Luminosity L l70 [h~ z L 



10" 



Fig. 6. The luminosity distribution of the FIRBACK 
170 /jm low-z galaxy sample. The dotted histogram is the 
luminosity distribution of the whole sample. It has an ex- 
cess at the highest luminosity bin. The dashed and dot- 
dashed histograms are those of the subsamples in South 
Marano field and ELAIS Nl field, respectively. 
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Fig. 8. Comparison of the 170-/im LF with those at 
other FIR wavelengths. Open triangles are our 170-//m 
L F (0 = 0.0), open diamonds represent the 100-/um LF 
of ISoifer fc Neugebauer Jl99lh. and open s quares repre- 
sent the 850-/im LF of iDunne et all <|2000lh The dashed 
curve is the analytic form of the 60-^tm LF of the 
IRAS PSCz entire sample, and the solid curve is the 
LF of the cool subsample of the IRAS PSCz galaxies 
l|Takeuchi. Yoshikawa. fc Ishiill2003[) . The 850-/im LF is 
horizontally shifted by a factor of 100 for display purposes. 



a LF at 100 /im from IRAS galaxy sample. Their 100- 
IMa LF has a bright end slightly steeper th an 60-/xm LF 
does, b ut the slope is still flatter than our LF. lDunne et alJ 
(2000) constructed a LF at 850 /mi based on 60-^m se- 
lected IRAS galaxy sample. Although the dynamic range 
is small, the over all shape of their 850 - urn LF is similar to 
our LF. Actually. IDunne et all (|2000h fitted their LF with 
the Schechter f unction which h as a very steep decline at 
the bright end <|Schechterl ll976). 6 

In summary the 170-/im LF has a shape similar to 
those of galaxy sample with cool dust emission. It is also 
interesting to note that the 170-/iin LF has the highest 
normalization among known FIR LFs. We summarize the 
comparison in Figure |H1 However, we must keep in mind 
the large uncertainty of the estimates, and further obser- 
vational exploration is definitely required. 

5.2. Evolution 

5.2.1. Pure luminosity evolution assumption 

Most of the galaxies in our redshift sample are at z < 0.3. 
Hence, we should estimate the evolution of the LF under 
a certain assumption. We adopt a pure luminosity evo- 
lution (PLE). Recent Spitzer observations indicated that 
the actual evolution of IR galaxies is described as a strong 
evolu tion in luminosity, with a slight evolution in densit y 
(e.g., IPerez-Conzalez et al .1 EooilLe Floc'h et ablEof)^. 
while their studies are based on the Spitzer 24 /jm band. 
Hence, the PLE is not a bad choice as a first approxima- 
tion. 

Adopting the PLE, the strength of the evolution can 
be estimated via a ra dial density distribution of galaxies 
l|Saunders et al.lll990j) . In the case of the PLE, the LF at 
z is expressed by the evolution strength f(z) as 



ML, z) = 4> 



L 

W) 



(11) 



where 4>o(L) is the local functional form of the LF. The 
PLE assumes that this form remains unchanged and only 
shifts along the luminosity axis. To define 0o(£i7o)> we 
examined the 170-/im LF for a sample at z < 0.15 (38 
galaxies). Though the uncertainty is large, we did not ob- 
serve a significant change of the LF in its shape, i.e., the 
PLE assumption might be approximately valid. Hence, we 
can use the shape of the LF of the sample with z < 0.3 
as 4>o(Lno). For simplicity, we consider a power-law form 
for f(z): 



f(z) = (l + zf 



(12) 



However, we sh ould keep in mind that the sample of 
IDunne et alJ i2000f) is selected at 60 /im. As these authors 
discussed, a significant fraction of galaxies with cold dust may 
be missed in the sample, and consequently, it might be possi- 
ble that their LF shape is not representative of cool/cold dust 
galaxies. 
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but this form is supported by recen t Spitzer observations 
at a wide redshift range of z = 0-1 ( Perez-Gonzalez et alJ 
l2005ULe Floc'h et alJl20ol . 



5.2.2. Estimation of evolution via radial density 

Assuming that the LF is separable for L and z as (j>{L, z) 
n(z)p(L), the likelihood is written as 



M(Q\{Li,Zi}i 



N 

n 

»=i 



4>{Li,Zi) 



d 2 V 
4>{Li,z) ,„ , dz 



(13) 



dfldz 



where N is the total sample size, z m ; n is the lowest red- 
shift which we used in the analysis, z maxi is the maxi- 
mum redshift to which ith galaxy can be detected, and 
(d 2 V/dildz)dz is again the differential comoving volume 
[Eq. (H3J]. 

By maximizing Equation (|13[) with respect to Q, we 
can have a maximum likelihood estimate. We performed 
this procedure with the nonparametric LF, because as 
we already mentioned, the analytic form underestimates 
the bright end of the LF, and it would lead to a seri- 
ous overestimation of the evolution strength. We obtained 
Q = g. The expected redshift distribution with this 

evolution is obtained by a similar manner to Equation © 
as 



dN_ 

dz 



n 



survey / S^O 
logL min (z) 



^170 



1 + 2 



d 2 V 
dQdz 



dlog£i 70 (;4) 



and shown in Figure[7| (the thin solid line). The agreement 
with the data and the expected value is much improved. 

Although the uncertainty is very large because of 
the limited sample size, this is similar to the value ob- 
tained by recent Spitzer 24-/im observations, Q ~ 4 
l|Perez-Gonzalez et al-lEooltlLe Floc'h et al.ll2005|) . These 
authors also found a weak evolution in the galaxy number 
density, but for the present sample, it is impossible to ex- 
plore this effect. If confirmed, the similarity between the 
strength of the galaxy evolution in the MIR (12-/um in the 
rest frame) and FIR will provide us an important clue to 
the physics of dusty star formation in galaxies. 

5.3. FIR luminosity density and obscured star 
formation density in the Local Universe 

5.3.1. The local 170-jum and total IR luminosity 
density 

The luminosity density in a cosmic volume provides vari- 
ous information of the energy distribution in the Universe. 
Especially, comparison between the radiative energy di- 
rectly emitted from stars and that re-emitted from dust is 
one of the key quantities to understand the fraction of hid- 
den star formation. In this subsection, we discuss the local 
luminosity density at 170 fim, pi(170^im) [L Q Mpc -3 ] and 
consider the integrated SED of the Local Universe. 
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Fig. 9. The MIR and FIR luminosity densities in the Local 
Universe. The filled square represent the 170-/im luminos- 
ity density, /5i(170/im), obtained by a direct integration 
of our LF. The vertical error bar shows the total uncer- 
tainty related to the LF. The horizontal bar represents 
the redshift range of our local sample (0 < z < 0.3). The 
filled triangle is the corrected pL(170^m) by assuming a 
pure luminosity evolution L(z) oc (1 + z) 5 . Open squares 
present the luminosity densities at 12, 25, 60, 100, and 
850 [im calculated based on the literatures (see text). 



In principle, it is straightforward to obtain pi(170/im) 
from the LF: we simply integrate the first-order mo- 
ment of the LF, Li7o0(-Li7o), over the whole possi- 
ble range of luminosity. Since the luminosity range of 
the 170-^tm LF is limited to 10 9 -10 12 L , we extrap- 
olated the faint end without observed data. We have 
done it by using the analytic form (Eq. EJ) with vari- 
ous a. However, as far as a < 2.0, the integration of 
ii7o0(ii7o) converges and the faint end does not con- 
tribute to the total integration significantly. For the in- 
frared _JlT{j_gahixies^ previous studies suggest a < 2.0 
(e.g., ISaunders et~aIlll99fA iTakeuchi. Yoshikawa. fc Ishii 
2003, and references therein), and the pL(170^m) is lit- 
tle affected by a. By the same reason, the lower and up- 
per bounds of the integration do not affect the result. We 
chose 10 7 /i _2 [L©] as the lowest luminosity, and use the 
highest luminosity bin as the upper bound. Thus, we ob- 
tainedp L (170Aim) = (1.1±^;|) x 10 8 /i[L Q Mpc~ 3 ]. The final 
uncertainty of yOi(170/xm) is dominated by the statistical 
uncertainty of the nonparametric LF in each bin. 

We, however, must recall that the present FIRBACK 
sample is not exactly 'local', i.e., it consists of galax- 
ies < z < 0.3. As seen in the previous subsection, it 
may be plausible that this result is affected by the strong 
galaxy evolution. If we adopt a PLE L(z) oc (1 + z) 5 , 
Pi(170^m) should be enhanced by a factor of 2.1 on av- 
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Table 1. The SED of the luminosity density in the Local 
Universe. 



A 


PL (A) 


Pl(A)/pl(170 p.m) 


H 


h [L Mp C - 3 ] 




12 


1.79 x 10' 


3.3 x 10" 1 


25 


1.42 x 10 7 


2.6 x 10 _1 


60 


4.08 x 10 7 


7.5 x 10" 1 


100 


8.25 x 10 7 


1.5 


170 


5.20 x 10 7 a 


1 


850 


9.45 x 10 4 b 


1.7 x 10" 3 



a The effect of the evolution is corrected assuming a pure 
luminosity evolution with L(z) oc (1 + z)® (Q = 5.0). 

b This value is calculated by integ rating a Schechter function 
presented bv lDunne etalJ«2000), over the luminosity range 
of L 850 = 1O 7 L -1O 13 L . 



erage in this redshift range. By correcting the evolution, 
we have p L (170^m) = 5.20 x 10 7 ft[£©Mpc _3 ]. 

We plot the result in Figure The filled square 
represents the 170-p.m luminosity density obtained by 
a direct integration of our LF. The horizontal bar in- 
dicates that the redshift range of our law-z sample is 
< z < 0.3. The filled triangle shows the evolution- 
corrected value. Then, we consider the SED of the lumi- 
nosity density in the Local Universe. We do similar ex- 
ercise s to calculate pj ( A) a t various IR wavele ngths. At 
MIR. lFang et alJ l)l998l) and lShupe et all l|l998|) provided 
the LFs at IRAS 12 and 25 mum bands, respectively. 
iFang et "ail l)l998[) tabulated their nonparametric LF at 
approximately the same luminosity range as we adopt 
in this work (1O 7 L Q -1O 12 L0). We simply summed it up 
with multiplying the lumin osity and obtained pl(12/iiti) = 
1.79 x 10 7 /i [L^Mpc'^. IShupe et all l)l998)) provided an 
analytic fit for their 25-/im LF. By their analytic function, 
we g ot p L (25pm) = 1.42 x 10 7 ft [L^Mpc " 1 ]. At 60 pm, 
from iTakeuchi. Yoshikawa. fc Ishiil ll2008h. we obtained 
p r ,(60 um) = 4.08xl0 7 ft [L^Mpc-^. ISoifer fc Neugebaueil 
(1991) presented the LFs at all the IRAS bands. Using 
their 100-pm nonparametric LF, we re-calculated the 
luminosity density to obtain pi(100/im) = 8.25 x 
10 7 ft fL^Mpc" 1 ]. Lastly, we used t he Schechter function 
fit provided bv iDunne et alJ l|2000l) to have pL(850/mi). 
Since their faint-end slope is very steep (a = 2.12), the 
integration is dependent on the adopted lowest luminos- 
ity in this case. We coherently integrated the Schechter 
function in the same range as the other bands. We found 
p L (850fim) = 1.00 x 10 5 /i [LqMpc" 1 ]. 

We plot these luminosity densities in Figure (open 
squares). The overall peak of the SED of the luminos- 
ity density seems to lie at A ^ 100 /urn. Further, if the 
suggested strong evolution is true, the local SED peak is 
restricted to be at 100 <; A ^ 170 /an. Even though the 
evolution effect significantly reduces the local pi(170/im) 
value, it still considerably contributes to the total FIR 
luminosity density. 



To have a crude estimate of the total IR luminosity 
density (we call it pi(dust)), we logarithmically inter- 
polate and extrapolate between the SED data points in 
units of [ergs -1 , Hz _1 Mpc~ 3 ], and integrate it over the 
range of 8-1000 p,m to match the conventiona l definition 
of th e total IR (TIR) luminosity (see, e.g., Dale et all 
EoOlllPale fc HeloulEool ITakeuchi et al.l l2005lh Tf we do 
not assume the evolution, we have pi(dust) = 1.4 x 
10 8 ft [L Q Mpc~ 3 ], a nd with the evolution, p r ,(dust) = 1-1 x 
10 8 ft [L Mpc- 3 ]. ITakeuchi. Buat. fc Burgarellal l|2005l) 
estimated the pl (dust) under the assumption of a con- 
stant ratio of p^(d ust)/pr,(60/im) = 2.5 (see, e.g., 
ITakeuchi etaP l2005|) . They found p L (dust) = 1.02 x 
10 8 /i [L Mpc -3 ]. Our pi (dust) values are in a very good 
agreement with that, especially for the case with evolu- 
tion. 

For the interpretation of the IR luminosity density, 
it should be worth mentioning that there is a possibly 
high contamination by active galactic nuclei (AGN) at 12 
and 24 pm. However, by integrating over the energy den- 
sity from these wevelengths, we find that the contribution 
from 12 and 24 p.m bands to pl (dust) is less than 10 %. 
Hence, if all the energy from 12 and 24 jxm were from the 
AGN, the effect of AGNs would not change the physical 
interpretation of pz,(dust) significantly. 

5.3.2. The obscured star formation density in the 
Local Universe 

The ratio between the energy from young stars directly 
observed at UV and that reprocessed by dust and observed 
at IR in the cosmic history has long been a matter of 
debate. Before closing the discussion, we consider the star 
formation rate density in the Local Universe obscured by 
dust. We use the value with evolutionary correction in the 
rest of this work. To get values without this evolutionary 
correction, we may simply substitute the former value for 
Pl (dust). 

For the conversion from pi(dust) to the cosmic star 
formation rate (SFR) densit y related to dust, psFit(dust), 
we can use several methods. LKennicutd lll998|) presented a 
famous formula between the SFR and the dust luminosity, 



SFR [MQyr- 1 ] = 1.72 x lO -10 L dust [L ] , 



(15) 



which is valid for starburst galaxies with a burst younger 
than 10 8 yr. Adopting this formula to pi(dust), we ob- 
tained p S FR.(dust) = 1.8 9 ft x IP" 2 [M ? yr^Mpc" 3 ]. 
However, as mentioned bv iKennicutd £l998) himself, this 
formula is not valid for more quiescent, normal galaxies. 
Since our SED of the Local Universe is similar to a kind 
of cool galaxies, we should caref ully treat the effect of the 
heatin g radiation from old stars. lHirashita. Buat. fc Inouel 
(2003J found that about 40 % of the dust heating in 
the nearby galaxies comes from stars older than 10 s yr. 
If we apply the correction of old stellar population 
to Equation (JTSJ, we obtain psFFt(dust) = 1.14 ft x 
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1(T 2 [M Q yr- x Mpc~ 3 ]. Belli (|2003l) also presented a simi- 
lar correction factor (0.32 ± 0.16 for galaxies with Ldust < 
10 n L Q and 0.09±0.05 for those with L dust > 10 11 J L Q ) for 
the contributio n of old stars. Based on a more theor etical 
point of view, iTakeuchi. Buat. fe Burgarellal l)2005|) also 
obtained an appropriate formula including the correction, 
which can be written as 

SFR [M Q yr" 1 ] = 1.07 x lO" 10 i dust [L Q ] . (16) 

Adopting Equation H16[l. we obtain psFR,(dust) = 
1.17 ft x 10~ 2 [M yr -1 Mpc~ 3 ], very close to the 
above. The obtained psFR(dust) is slightly larger 
than the local SFR density estimated from direct 
FUV radiation (without dust attenuation correc- 
tion), gSFR (FUV) = 8.3 x 10^[M Q yr^Mpc" 3 ] 
(ITakeuchi. Buat. fc Burgarellal 1200.4 see also 
ISchiminovich et all 1200.4) . Hence. 58.5 % of the 
star formation is obscured by dust in the Local 
Uni verse. This is in very good agree ment with that 
of ITakeuchi. Buat. fc Burgarellal (|200.4) . but since we 
reached this conclusion from the measured dust SED 
of the Local Universe, we could put a firmer basis on 
their conclusion by this work. These results may be the 
first direct estimate of the dust luminosity in the Local 
Universe, and should be tested by forthcoming large area 
survey in the FIR by e.g., ASTRO-F. 7 



6. Conclusion 

We analyzed the FIRBACK 170 pm galaxy sample to ob- 
tain the local luminosity function (LF) of galaxies. We 
constructed a flux-limited sample with Sno > 0.195 Jy 
and z < 0.3 from the survey, which consists of 55 galax- 
ies. 

The overall shape of the 170-^m LF is 
quite different from that o f the total 60-pm LF 
^Takeuchi. Yoshikawa. fc Ishiil |2003|) : the bright end 
of the LF declines more steeply than that of the 60-/im 
LF. This behavior is quantitatively similar to the LF 
of the cool subsample of the IRAS PSCz galaxies. The 
brig ht end is al s o sim ilar to that of the submillimeter LF 
of lDunne et"aD l|2000l) . 

We also estimated the strength of the evolution of the 
LF by assuming the pure luminosity evolution L(z) oc 
(1 + z) Q . We obtained Q = 5.0^ which is simi- 
lar to the value obtained by recent Spitzer observa - 
tions l)Perez-Gonzalez et al J 120051 iLe Floc'h et al.ll2005|) . 
in spite of the limited sample size. 

Then, integrating over the 170-pra LF, we ob- 
tained the local luminosity density at 170 pm, 
Pi(170^m) [L Mpc~ 3 ]. If we assume the above strong 
luminosity evolution L(z) oc (1 + z) 5 , the value is 5.2 x 
10 7 ft [L Mpc~ 3 ], which is a considerable contribution to 
the local FIR luminosity density. 

By summing up the other MIR/FIR data, we ob- 
tained the total dust luminosity density in the Local 



7 URL: |http : //www . ir . isas . ac . jp/ASTRO-F/index-e .html 



Universe. We obtained p L (dust) = 1.1 x 10 8 ft [i Mpc 3 ] 
without the evolution correction, and p^(dust) = 1.1 x 
10 8 ft [£ Mpc~ 3 ] with correction. 

Lastly, based on pl (dust), we estimated the cos- 
mic star formation rate (SFR) density hidden by dust 
in the Local Universe. We took into account the dust 
emission heated by old stellar population, and ob- 
tained p SFR (dust) ~ 1.1-1.2 ft x 10" 2 [Moyr^Mpc -3 ]. 
Comparin g with the SFR density estimated for m FUV ob- 
servation (Tak euchi. Buat. fc Burgarellall2005|) . we found 
that 58.5 % of the star formation is obscured by dust in 
the Local Universe. 

It will be important to examine our local LF by a large 
area survey of the Local Universe. The ASTRO-F project 
promises to provide a local large sample of FIR galaxies 
at ~ 50-170 pm. For the evolutionary status of the FIR 
galaxies, the FIR data of Spitzer will be very important 
to examine the present result. 
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Fig. B.l. Comparison of the luminosity functions of our 
ISO 170 /mi galaxy sample by 1/Vmax and C~ estimators. 
Thick and thin symbols represent the C~ nonparametric 
LFs, respectively. The signification's of the symbols are 
the same as Figure 5 The symbols for 1/U max estimates 
are shifted with 0.05 dex for th purpose of visual clarity. 
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Appendix A: Tables of the nonparametric 
luminosity function 

In this section, we present numerical tables of the nonpara- 
metric luminosity function (LF) of the FIRBACK 170 fim 
galaxy sample obtained with the C~-method. We only 
tabulate the LF directly estimated from the original data, 
and the flux measurement errors are not included in these 
LFs. 

Appendix B: Comparison between the l/V mSLX - 
and C -luminosity functions 

In this section, we examine the problem of the classical 
1/Vmax estimator (Schmidt 1968; Eales 1993). Comparison 
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Table A.l. The 170-fim luminosity function with j3 = 1.0, estimated with Lynden-Bell's C~ method. 



logLi7o a 


(j)(Ll7o) 

h 3 [Mpc^dex" 1 ] 


0(ii7o) upper 
h 3 [Mpc-Mex" 1 ] 


J / r \ lower 

h 3 [Mpc^dex" 1 ] 


9.15 


4.89 x 10 -2 


6.20 x 10 -2 


0.00 x 10 +u 


9.45 


3.10 x 10 -2 


4.96 x 10" 2 


1.24 x 10 -2 


9.75 


1.73 x 10" 2 


2.63 x 10" 2 


7.70 x 10" 3 


10.05 


3.81 x 10~ 3 


5.99 x 10~ 3 


1.43 x 10 -3 


10.35 


1.81 x 10~ 4 


3.00 x 10~ 4 


5.95 x 10~ 5 


10.65 


2.29 x 10 -5 


3.85 x 10~ 5 


5.54 x 10~ 6 


10.95 


6.87 x 10" 6 


1.17 x 10" 5 


0.00 x 10 +0 



a Units of Li7o is h 2 [ L@], and the bin width is 0.3 dex. Tabulated log luminosities represent the bin center. 
Table A.2. Same as Table lAdl but for (3 = 0.0. 



log L 170 


4>{Lno) 


0(Li7o) upper 


±(r \ lower 
0(^170 J 




h 3 [Mp C - 3 dex- 1 ] 


h 3 [Mpc- 3 dex- 1 ] 


h 3 [Mpc- 3 dex- 1 ] 


9.15 


3.33 x 10" 2 


4.65 x 10~ 2 


0.00 x 10 +u 


9.45 


5.12 x 10" 2 


7.17 x 10~ 2 


3.10 x 10" 2 


9.75 


1.69 x 10~ 2 


2.56 x 10" 2 


7.42 x 10~ 3 


10.05 


4.29 x 10~ 3 


6.80 x 10" 3 


1.69 x 10~ 3 


10.35 


1.89 x 10~ 4 


3.16 x 10~ 4 


5.95 x 10~ 5 


10.65 


4.04 x 10 -5 


7.00 x 10" 5 


7.20 x 10~ 6 


10.95 


1.20 x 10~ 5 


2.01 x 10 -5 


0.00 x 10 +0 



Table A.3. Same as Table fO but for /3 = -0.5. 



logLl70 


4>{L 170 ) 
h 3 [Mpc- 3 dex- 1 ] 


0(Li 7O ) upper 
h 3 [Mpc- 3 dex- 1 ] 


i ( r \ lower 
</>(-Ll70j 

h 3 [Mpc- 3 dex- 1 ] 


9.15 


3.37 x lGT* 


4.65 x 10~' z 


0.00 x 10 +u 


9.45 


5.14 x 10~ 2 


7.21 x 10~ 2 


3.10 x 10 -2 


9.75 


1.67 x 10 -2 


2.54 x 10~ 2 


7.45 x 10~ 3 


10.05 


4.12 x 10 -3 


6.57 x 10" 3 


1.53 x 10~ 3 


10.35 


1.56 x 10 -4 


2.63 x 10 -4 


4.51 x 10~ 5 


10.65 


5.45 x 10 -5 


9.37 x 10" 5 


1.11 x 10~ 5 


10.95 


1.40 x 10~ 5 


2.40 x 10 -6 


0.00 x 10 +0 



of the luminosity functions of our ISO 170 /im galaxy sam- 
ple by 1/Vmax and C~ estimators is shown in Figure lB~TI 
In Figure iRll the l/V mayi LFs are shifted to the left with 
0.05 dex to make them easy to see, but the actual bin 
centers are exactly the same as those of C~ LFs. It is im- 
pressive that there is a large difference between the 1 /V mayi 
and C~ LFs. At the fainter side, these two LFs are con- 
sistent with one another within the error bars, though 
the faint end of the 1/Vm ax LFs tend to be slightly un- 
derestimated. In contrast, the bright end is completely 
different: 1/Vinax method gives much flatter LFs than C~ 
method. We should recall that the parametric method and 
the C~ method are both insensitive to density fluctua- 
tion, while 1/Vmax method is unbiased only for the spa- 
tially homogeneous sample, as extensively examined by 
Takeuchi et al. (2000). Then, the most plausible explana- 
tion of this discrepancy may be due to the existence of 
a density enhancement, corresponding to the luminosity 
L 170 ;> 3 x 10 10 Lq. 

To understand more clearly, we recall the redshift dis- 
tribution of the present sample (Figure 7). The distribu- 



tion of ELAIS Nl galaxies is not very far from that ex- 
pected by the LF, while that of South Marano galaxies 
are very different. Reflecting the luminosity distribution 
of this field, which is heavily inclined to luminous galaxies, 
the distribution has no apparent peak, but has a widely 
spread shape toward redshifts up to ~ 0.3. This heavy 
high-z tail is superposed to the tail of the ELAIS Nl field 
makes a significant bump at z ~ 0.2. Since the 1/Knax 
method assumes a spatially homogeneous source distri- 
bution, this excess is too much exaggerated through the 
number density estimation, and results in the overestima- 
tion of the corresponding luminosity bins. 



